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Using first-principles methods based on density-functional theory we investigate the spin relaxation in W(OOl) 
ultrathin films. Within the framework of the Elliott- Yafet theory we calculate the spin mixing of the Bloch states 
and we explicitly consider spin-flip scattering off self-adatoms. We find an oscillatory behavior of the spin- 
mixing parameter and relaxation rate as a function of the film thickness, which we trace back to surface-state 
properties. We also analyze the Rashba effect experienced by the surface states and discuss its influence on the 
spin relaxation. Finally we calculate the anisotropy of the spin-relaxation rate with respect to the polarization 
direction of the excited spin population relative to the crystallographic axes of the film. We find that the spin- 
relaxation rate can increase by as much as 47% when the spin polarization is directed out of plane, compared to 
the case when it is in plane. Our calculations are based on the multiple- scattering formalism of the Korringa- 
Kohn-Ro stoker Green-function method. 
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I. INTRODUCTION 

In the process of spin relaxation an excited electron spin 
population returns to the state of equilibrium that, in non- 
magnetic materials, corresponds to zero spin polarization. De- 
spite the fact that the fundamental mechanisms contributing 
to spin relaxation have been investigated since a long time in 
various systems, the phenomenon still attracts attention, ow- 
ing to its importance in spintronics applications,^ for example 
in giant magnetoresistance or, lately, in the spin Hall effect as 
well as the inverse spin Hall effect that is used to probe spin 
currents.^ We mention these examples among a variety of ap- 
plications in order to stress the practical importance of spin 
relaxation in thin metallic films, which is part of the motiva- 
tion for the present work, as a source of loss of spin-mediated 
information. 

There are various mechanisms that can contribute to spin 
relaxation^ ^ and in metallic systems they are mostly related 
to Fermi- surface properties. Although it is clear that in a non- 
magnetic metal or metallic film the spin relaxation can be at- 
tributed to spin-flip scattering in the presence of the spin-orbit 
coupling (SOC), many parameters come into play, and one ex- 
pects the spin-relaxation rate to depend strongly on the film's 
crystallographic orientation, its thickness and details of the 
Fermi surface. Of particular importance here^ can also be the 
Rashba states created as a result of spin -orbit interaction act- 
ing on the surface bands in the film.tSE^ 

In this work we present a study of free-standing W(OOl) 
films with space-inversion symmetry where we consider that 
the spin-flip scattering is induced by W adatom impuri- 
ties. The motivation for restricting the investigation to free- 
standing films is that we consider them to be generic proto- 
types for films in layered structures if the material "sandwich- 
ing" the film — the surrounding matrix — is insulating. There 
is no question that contact at the film surfaces will produce ef- 
fects that depend on the surrounding matrix, however, we are 
searching here for physical mechanisms that can in principle 



still be present in the contact case, even if the details of the 
band structure are different. In the present analysis we will 
face e.g. oscillatory effects arising from the thickness depen- 
dent interaction of surface states that can be replaced by inter- 
face states if the films are sandwiched; also, the low dimen- 
sionality entails an anisotropy of spin relaxation that should 
be a quite general effect irrespective of the film contact. Due 
to the pronounced manifestation of such effects, bcc W(OOl) 
films are chosen among bd transition-metal thin films for a 
deeper analysis in the present work. 

In the presence of structural inversion symmetry, which 
is the case here in bcc(OOl) films, the Elliott- Yafet 
mechanism^ ^ plays the most important role for spin relax- 
ation. Within this mechanism, the relaxation is realized via 
spin-orbit mediated spin-flip scattering off impurities at low 
temperatures and additionally off phonons at higher tempera- 
tures. According to Elliott,^ in a system with time-reversaP^ 
and space-inversion symmetry there are two degenerate Bloch 
states at each k-point, which can be written as superpositions 
of up I t) and down | \) spinors: 

*+(r) = (ak(r)|t)+6k(r)U»e''^-'- 
*k (r) = (alk(r)l \) " ^IkWI t)) e*-. 

Due to the degeneracy there is an arbitrariness in the selec- 
tion of ak(r) and 6k (r), as any superposition of the states in 
Eq. ([T]) is also an eigenstate of the Hamiltonian. In practice 
the arbitrariness is usuall>^ lifted by the demand that the 
spin expectation value of should be niaximal in the z di- 
rection, or actually in any chosen directioiP321 s (that we call 
the spin-quantization axis, SQA) experimentally defined by 
the polarization direction of the excited spin population. The 
spin polarization vector corresponding to can be cal- 
culated via 

= ^(*^kl*^)> i = x,y,z, (2) 

where cr is the vector of Pauli matrices (atomic units with h = 
1 are implied). The two degenerate wavefunctions at each 
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k-point are orthogonal to each other and have opposite spin 
expectation values: 
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It is straightforward to show the relation of the spin projection 
along the SQA, 5 • Sj, to ak(r) and 6k (r): 
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where we have also introduced the integrals and 6^. Taking 
an average over the Fermi surface (FS) (actually Fermi lines in 
the present two-dimensional case), the Elliott- Yafet parameter 
(EYP) is obtained as 

= (^k)pq = ,r^\,r I dk^—, (6) 
\ ^/FS n(^F)^Bz7FS ^F(k)' 

where VF(k) is the Fermi velocity, n{EF) is the density of 
states at the Fermi level and Vbz is the Brillouin zone 
volume. 

Momentum scattering events couple the spin-up and down 
components of the wavefunctions at different momenta, al- 
lowing for transitions between ^ J and and giving rise 
to spin relaxation. In the Elliott approximation^ the spin- 
flip probability rate of state ^J^, , is proportional to the 
momentum-dependent spin-mixing parameter 6^. As a re- 
sult, after taking an average over the Fermi surface, the ra- 
tio between the spin-relaxation rate 1/Ti and momentum- 
relaxation rate 1/Tp, T^^/T~^, is proportional to 6^. De- 
parting from the Elliott approximation, in which P is as- 
sumed small, one has to take into account the form of the 
scattering potential of the impurity (W adatom in this 
case) and calculate the scattering T-matrix, T{E) = AV(1 — 
G{E)AV)~'^, where G{E) is the Green function of the un- 
perturbed film. Then, the spin-conserving and spin-flip prob- 
ability rates are given by the squared matrix elements = 
27r|(*+|r(^F)|*+)pandPi7 = 27rK*+|T(^F)|*- )p, 
respectively. Fermi-surface integrals of the scattering proba- 
bility yield the average momentum- and spin-relaxation rate 
per k-point 
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with c the concentration of impurities. The factor 2 in the def- 
inition of Ti with respect to the spin-flip time Tgf has a histor- 
ical origin as Ti was derived from the full linewidth at half- 
amplitude of conduction elecron spin resonance spectra.^^I^ 



In the literature, Ti is usually mentioned as the measured 
quantity, but from the point of view of scattering theory it is 
more natural to use Tgf . 

It is clear that the spin-mixing parameter 6^ reflects the host 
contribution to the spin relaxation while the spin-relaxation 
time Ti contains also the scattering contribution of the impu- 
rity which of course depends on the type of impurity as well 
as impurity concentration. In the present work, we will first 
discuss the spin-mixing parameter to understand general fea- 
tures of spin relaxation in W(OOl) thin films. After that, by 
introducing a W adatom on one film surface, we will inves- 
tigate the spin-relaxation rate T^^ quantitatively. The spin- 
quantization axis s is at first chosen perpendicular to the sur- 
face of the films, but later we also examine a variation of the 
quantization axis revealing anisotropic effects in spin relax- 
ation. The variation of the quantization axis corresponds to 
an experimental situation where the spin polarization of the 
injected spin population (or the magnetic field direction in an 
electron spin resonance experiment) is changed. 

We find that the EYP acquires very large values owing to 
the Rashba effect at the surface and also exhibits an even- 
odd oscillation with the thickness of the films following the 
behavior of the electronic structure of the surface bands; the 
same is true for the spin-relaxation rate due to scattering off 
W adatoms. These effects are the subjects of Sec.[ll|and III 



We then present in Sec. IV the anisotropy of the EYP and of 
the spin-relaxation rate with respect to the choice of the SQA. 
Finally in Sec.|V|we argue that, as far as the spin relaxation in 
the surface states is concerned, the mechanism discussed and 
calculated here should be dominant over, e.g., the Dyakonov- 
Perel mechanism, prominent in semiconductors or semicon- 
ductor heterostructures where Rashba-type of splitting is also 
present"^ ^ but has a much lower magnitude. 

Our investigation is based on density-functional calcula- 
tions within the local density approximation.^^ We employ the 
full-potential Korringa-Kohn-Rostoker (KKR) Green function 
method^ with exact treatment of the atomic cell shapes. 
After a self-consistent full-potential calculation performed 
within the scalar-relativistic approximation, spin-orbit cou- 
pling is added when calculating the Fermi surface properties 
and the scattering-matrix elements. The formalism for the cal- 
culation is given in detail in Ref. |9|. A similar formalisnP^ 
has been applied before in Refs. 1221231 where the results 
compared well with experiment. An angular momentum cut- 
off of /max = 3 is taken. We use the experimental lattice 
constant of bcc tungsten of a = 3.165 A. For an interpreta- 
tion of our results and a separation of causes we perform in 
some cases numerical experiments by switching on and off 
the spin-orbit coupling in the system. The W adatom im- 
purity is treated using the Jiilich KKR impurity-embedding 
code (KKRimp) which enables us to treat with the charge and 
spin-density self-consistently including the perturbation of the 
nearest neighbors of the impurity. Lattice relaxations are not 
taken into account. 

It should be mentioned that W(OOl) films have been a sub- 
ject of intensive theoretical and experimental reasearch for 
a long time, cf. Refs. lf24K26l and citations therein. It is 
known that at low temperatures the W(OOl) surface under- 
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FIG. 1: (Color online) The atom-resolved density of states (DOS), 
calculated ignoring the SOC, of the surface- and central layer of a 7- 
layer thick W(OOl) film in comparison to the DOS of bulk W shown 
as a grey- shaded area. 



goes a c(2 x 2) reconstructio n, wh ile the ideal structure is re- 
stored at higher temperatures P^^^ Additionally, the bcc struc- 
ture of W(OOl) films has been extensively studied owing to its 
multiple applications as a substrate used for deposition P^^^ 
Recently, because of the strong spin-orbit interaction W(OOl) 
is used as subtrate for ultrathin magnetic films generating a 
strong Dzyaloshinki-Moriya interaction introducing complex 
magnetic structures. However, most important for our study 
are the surface states. Together with Cu(lll), W(OOl) is one 
of the first syste ms for which surface states were predicted to 
exist theoretically.l3El] 
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FIG. 2: (Color online) Fermi surfaces calculated including SOC and 
distribution of the spin-mixing parameter bl- (color of the Fermi sur- 
face points) for W(OOl) films of various thickness. The Fermi sur- 
faces are grouped into those for films with even number of layers 
(upper group) and odd number of layers (lower group). The sur- 
face states SI and S2 are marked with arrows. Small circles indicate 
where S2 merges into the bulk-like states. 



II. EVEN-ODD EFFECT IN SPIN-MIXING PARAMETER 
AND SPIN-RELAXATION RATE 

We first present our results on an unexpected even-odd os- 
cillatory variation that we observed in the Elliott- Yafet param- 
eter of the film and the spin-relaxation rate due to a W adatom 
with respect to the number of layers in the W(OOl) film. 

To start the discussion, we take a look at the comparison 
of the density of states (DOS) of W bulk and a typical rep- 
resentative of the films studied in the following — a 7-layer 
W(OOl) film — shown in Fig.[T] In this particular calculation 
the spin-orbit coupling was not included. Characteristic of a 
transition metal with bcc structure, the bulk DOS shows a bi- 
modular behavior with a dip in the middle separating bonding 
from antibonding states, which for W lies at the Fermi en- 
ergy. This structure of the DOS is also clearly preserved in 
the center of the thin W(OOl) film. In contrast, the local DOS 
in the surface atomic layer of the film displays a peak at E^, 
which can be attributed to the presence of surface states. Ob- 
viously, this should result in a significant surface contribution 
to the Fermi- surface-dependent quantities and in particular the 
ElHott-Yafet parameter. The DOS of W(OOl) films of differ- 



ent thickness, considered later, are very similar to those in the 
case of a 7-layer film. 

Next, we include the spin-orbit coupling, choose the spin- 
quantization axis perpendicular to the film (s\\ [001]) and cal- 
culate the Fermi surfaces and the distribution of the spin- 
mixing parameter 6^ for W(OOl) films of varying thickness. 
The results are presented in Fig. [2] The most important fea- 
ture for our study are the surface states. These are indicated 
as "SI" and "S2," and pointed at by arrows in the vicinity of 
the M point on the Fermi surface. SI and S2 are present at all 
film thicknesses and we will discuss them in more detail be- 
low. Here we merely note for clarity that the S 2-line merges 
into bulk-like states at its two ends (indicated by two small 
circles in Fig.|2]), while SI does not; also, that SI appears as 
a double line due to SOC-induced Rashba splitting and due 
to the interaction between the two film surfaces, while S2 ap- 
pears as a single line because its partner is higher than in 
energy. 

The importance of the surface states for spin relaxation 
can be ascertained by looking at the Fermi- surface-integrated 
Elliott- Yafet parameter 6^ [Eq. ^] as a function of film thick- 
ness, additionally decomposed into surface and bulk contri- 
butions and presented in the top panel of Fig. |3] The latter 
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FIG. 3: (Color online) Top: The dashed lines show the thickness de- 
pendence of the Elliott- Yafet parameter in W(OOl) thin films with 
the SQA perpendicular to the film. The total EYP (squares) is de- 
composed into contributions from the surface states SI and S2 (tri- 
angles) and from the bulk states (diamonds) of the film. The full 
line (circles) shows the ratio between spin-flip rate and momentum- 
relaxation rate, T^'^ /T^-^ for the case of self-adatom scattering [W 
adatom on W(OOl)]. The even-odd effect is evident in the ratio fol- 
lowing very closely the oscillation of the Elliott- Yafet parameter. 
Bottom: SpHtting of the state S2' (without SOC) due to the inter- 
action between the two surface states of opposite surfaces (see Fig.H] 
for the band structure and for the position where the splitting of S2 
is calculated). The oscillatory behavior of the splitting as function of 
film thickness correlates with the behavior of the EYP and spin-flip 
rate. The lines are guides to the eye. 



decomposition was performed by integrating in Eq. ^ over 
only the surface-state part or only the bulk-state part of the 
Fermi surface, while keeping the denominator n(£^F) fixed at 
the value of the total density of states. It is striking that for 
all considered film thicknesses the overall contribution of the 
surface states to the EYP is comparable, or even larger than, 
the bulk contribution. It is well known that the amplitude of 
the surface states drops exponentially as a function of distance 
from the surface. Therefore, in the limit of large film thick- 
ness, the relative contribution of the surface states decreases 
and the bulk contribution becomes prominent. However, in the 
ultrathin films studied here (maximum 12 layers thickness), 
the surface states are more like quantum-well states and do 
not fully decay in the center of the film. Additionally, there 
are two more reasons why the bulk limit is not reached fast. 
First, the density of bulk-like states at is at a minimum as 
we saw in Fig. [l] thus the bulk contribution to 6^ sets in only 
slowly. Second, the Fermi surface of bulk W is rather complex 
with the value of 6^ strongly varying over it (see e.g. Ref. 14 ), 
thus many film layers are needed to achieve the equivalent of 
a good resolution in the [001] direction of the bulk Brillouin 
zone that would yield the bulk limit 6^ = 0.065.^^ Surpris- 
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FIG. 4: (Color online) The band structure of 6- and 7-layer W(OOl) 
films with (filled blue squares) and without (red open circles) spin- 
orbit coupling. Arrows mark the surface states SI and S2 correspond- 
ing to the calculation with SOC and Sl^ and S2^ corresponding to 
the calculation without SOC. An even-odd effect in the thickness- 
dependent splitting of S2^ is responsible for the even-odd effect in 
the EYP 



ingly, as we observe from Fig. [3] the contribution of the sur- 
face state 82 and correspondingly the total EYP displays a 
pronounced even-odd oscillation as a function of film thick- 
ness up to 10 layers. 

If we examine the distribution of on the Fermi surface 
for several film thicknesses, presented in Fig. [2j it becomes 
apparent that 6^ exhibits large variations in magnitude, both 
as a function of the position at the Fermi surface as well as the 
film thickness. We also clearly observe that among all states at 
the Fermi surface the largest spin-mixing occurs in the surface 
states SI and S2. 

Let us try to understand the oscillation with film thickness 
starting from the band structure and symmetry properties of 
the surface states. We consider the surface states in the ab- 
sence of spin-orbit coupling because it is simpler to apply 
symmetry arguments in this case. We name the states SI' and 
S2' where the primed indices are used to distinguish the case 
without SOC; SI' and S2' are mixed with each other under the 
action of the SOC Hamiltonian to produce SI and S2. 
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We consider the band structure of SI' and S2' along the FM 
diagonal of the Brillouin zone (red open circles in Fig.|4]). SI' 
and S2' comprise mainly d orbitals. We set a coordinate sys- 
tem with X and y along the [100] and [010] directions in the 
surface and z normal to the surface. We distinguish two types 
of d orbitals with respect to their reflection properties about 
the diagonal [110] (FM): even [^{d^z + dyz)^ dxy^ and dz^] 
and odd [^{dxz — dyz) and dx2_y2]. Staying on the high- 
symmetry line k G FM, the Bloch states ^k(r) derived from 
odd orbitals show nodes along the [110] diagonal, which has 
a consequence of a high in-plane kinetic energy, not leaving 
enough energy for penetration into the bulk region. Therefore, 
the odd-(i-orbital surface states of the opposite surfaces cou- 
ple very weakly to each other and show an almost vanishing 
splitting already at small film thicknesses. On the other hand, 
the even states, relieved from this nodal structure, have less 
in-plane kinetic energy and thus enough energy to penetrate 
into the bulk region and hybridize with their likes of the oppo- 
site surface; then the resulting hybrids show a splitting even 
at larger thicknesses. This parity-dependence of the surface 
states was observed in the past by Mattheiss and Hamann.^^^ 
The consequence of this can be seen e.g. in Fig.[4]for the bands 
without spin-orbit coupling along FM where the S2' bands 
show a large splitting whereas the bands SI' show an almost 
complete degeneracy. 

The even-odd effect of the EYP of state S2 can now be 
traced back to surface state S2'. The aforementioned splitting 
of S2' due to the interaction between the opposite surfaces ex- 
hibits an even-odd behavior: for an even number of layers the 
splitting is large, for odd it is considerably smaller due to an 
oscillation in the coupling. This can be seen in Fig. [3] (bottom 
panel) where the splitting of S2' has been calculated at the 
position indicated by arrows in Fig. [4] One immediately rec- 
ognizes the striking correlation between the oscillations in the 
splitting of S2' and in the EYP shown in the same figure. We 
should comment that the origin of the even-odd effect in the 
coupling could not be fully explained by the symmetry prop- 
erties of the wavefunctions. However, it is clear that owing to 
the different inversion symmetry center in even- and odd-layer 
films, the overlap between the opposite-side surface states is 
different which contributes to the even-odd effect. 

Including the SOC, the even-odd dependence of the en- 
ergy splitting will have a profound effect on the spin- 
mixing parameter, since the energy distance between states 
affects crucially the value of 6k in Eq. ([T]). To justify 
this statement we remind the reader of Elliott's perturba- 
tive expansion Pl for the n-th band a summation over all 
other bands n' ^ n should be performed, reading h'^^ = 

E„, |(*°kl^(i5)^^|*o,k>|V(^°k - K'^f in first order 
perturbation theory, where the superscript refers to the wave- 
functions and eigenstates without SOC, ^ is the SOC strength 
and {LSY^ is the spin-flip part of SOC operator. Since the 
energy difference of the states appears in the denominator (in 
all orders of the perturbation expansion), changing the energy 
splitting will strongly affect the value of 6^. Thus emerges the 
correlation between the splitting of S2' without SOC and the 
EYP with SOC in Fig.jS] 



Now we investigate the spin relaxation due to the W adatom 
impurities located on one film surface. In the upper panel of 
Fig. |3] we also show the ratio between the spin-flip rate and 
momentum-relaxation rate, T~^^ /T~^, as a function of the 
film thickness. We see that the magnitude of the ratio is of 
the same order as the Elliott- Yafet parameter and that the os- 
cillatory behavior of the two quantities is clearly correlated, 
even though Elliott's approximate relation T~^^ /T~^ ^ 46^ 
does not hold. The deviation from Elliott's approximation is 
not surprising since it holds under the assumptions that iP' is 
small and that the scattering is weak enough to be described 
within first-order perturbation theory. In any case the high val- 
ues of Tg7^/T~^, of the order of 0.3-0.4, show that approxi- 
mately every third scattering event includes a spin-flip. This 
high spin-flip rate is certainly related with the fact that a W 
adatom introduces an additional contribution to spin-flip scat- 
tering via its internal spin-orbit coupling. It can be remarked 
that the ratio between spin-flip rate and momentum-relaxation 
rate does not depend on the impurity concentration, but the 
values of the two quantities seperately do. In order to have 
a quantitative estimate, the impurity concentration c in Eq. [7] 
and [5] is set to be 1%. This gives us, for example, a value 
of the spin-relaxation rate of 1 1.63 ps~^/at% and momentum- 
relaxation rate of 13.69 ps~^/at% for a 10-layer W(OOl) thin 
film by scattering off W adatoms. We will return to the spin 
relaxation rate in SecHVl 



III. RASHBA-CHARACTER OF SURFACE STATES 

This section is dedicated to the investigation of the Rashba 
effect at the W(OOl) surface. The essence of the Rashba ef- 
fect at surfaces lies in a crystal-momentum dependent split- 
ting of surface states due to the presence of spin-orbit cou- 
pling together with the structural inversion asymmetry caused 
by the surf ace. ^^^^ As was first orbserved experimentally in 
Au(lll) by LaShell et al}^ and then theoreticaly explained 
by Henk and co-workers, the surface states of Au(lll) dis- 
play a SOC-derived dispersion relation very similar to that of 
the two-dimensional electron gas (2DEG) described by a two- 
band Rashba model. One signature of the Rashba effect is a 
Fermi surface consisting of two concentric rings with the spin 
polarization perpendicular to the group velocity but pointing 
in opposite directions at each of the rings. We will refer to 
the resulting k-space spin texture as spin-polarization field in 
the following (which is closely related to the spin-orbit field 
defined e.g. in Ref. 5 ). 

The Rashba splitting results in a degeneracy lifting where 
the pairing implied by Eq. ([T]) does not hold any more. Since 
the effect originates in the broken inversion symmetry at the 
surface, it can be captured computationally by treating half- 
infinite systems or by breaking the symmetry between two 
surfaces of a finite symmetric slab via deposition on a sub- 
strate or adsorption of a monolayer of a different material on 
one surface of the slab, e.g. a Cu monolayer on an Au filnP^ 
or H adsorption on W(llO).^^ Here, we follow an alternative 
approach, keeping the inversion symmetry but projecting out 
the appropriate state from the doubly degenerate subspace de- 
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FIG. 5: (Color online) Top: Schematic diagram showing the energy 
levels of the Rashba states at the film surfaces and their hybridization 
when the film thickness is finite and the states of opposite surfaces 
interact. The full lines represent the energy levels before and after 
hybridization and the dashed lines show which level pairs are hy- 
bridized. Bottom: The spin polarization of the surface states around 
the M-point of an 1 1 -layer W(OOl) film. The direction of the spin po- 
larization is given by arrows, while its magnitude is given by means 
of a color code of the arrows. Surface states SI and S2 are marked. 



fined in Eq. 

We introduce an intuitive view of our approach by means 
of a schematic energy level diagram in Fig. [5] (top). We first 
imagine the film in the limit of infinite thickness, with the 
states of each surface not interacting with the opposite sur- 
face. Then, a Rashba splitting occurs due to SOC. Fixing 
a k-point on the surface band, the splitting results in spin- 
polarized levels |L X) ^nd |L |) on the left surface (indi- 
cated by L) and similarly |R X) ^nd |R |) on the right sur- 
face. The Rashba splitting A is the same in the two surfaces 
but the order with respect to energy in which the levels of 
spin polarization t and | occur is reversed in the two surfaces 
because of the opposite direction of the surface normal, ac- 
cording to standard Rashba theory. Reducing the film thick- 
ness, the left and right states interact via a hopping t (let us 
assume t <C A) and hybrid states are formed with bonding 
(B) and antibonding (A) nature: |B t) - |L t) + t), 



|A t) ^ -i|L t) + |R t). Thus, |B t) is localized more 
on the left surface and | A X) more on the right surface. Anal- 
ogously, hybrids |B |) and |A |) are formed. A calculation 
of the film band structure finds linear combinations of the de- 
generate bonding a|B t) + /3|B |) as well as antibonding 
a| A t) + /3| A I) hybrids; such linear combinations are basi- 
cally the degenerate wavef unctions described 



by Eq. (m) 



in the particular case of surface states. Returning to the 
notation of Eq. ([T]), it has then to be decided which particular 
linear combination is of interest for the physics of the prob- 
lem at hand. For instance, for the calculation of the EYP in 
the previous section the linear combinations were chosen so 
that the spin expectation value along the z direction was max- 
imized. Here, on the other hand, we want to choose constants 
Qfk and (not depending on r) in such a way that the com- 
bination Q^k^k + /^k^k resembles as much as possible the 
Rashba states of the infinite-thickness film; i.e., in a way that 
the wavefunctions |B t), |B |), |A t), and |A |) of Fig.[5]are 
retrieved. For this reason we pick one surface, say the left, and 
we find the linear combination that maximizes the spin expec- 
tation value within the particular surface atomic layer and call 
the resulting wavefunction ^j^^^'^ = ak^k + /^k^k • 
thus define the polarization in the surface layer as 



k/surf 



max.L 



(r) d^r, 



surf .layer 



(11) 

G {x^y^z}, and demand that o^k and /3k are such that 



5k/surf I 



k / surf 



(SI)' 



surf 



(^k)surf is maximized. 



Then we observe that ^|^^^'^ has both the spin and charge 
density more localized on the left surface (numerical result not 
shown here explicitly). Since the film potential in our calcu- 

= a* 



lation is still inversion- symmetric, ^|^' 
has an orthogonal degenerate partner, ^|^^^'^ 
/3* k^k ' which is more localized on the opposite surface and 
has the opposite spin expectation value. We say "more lo- 
calized" and not "completely localized" because ^j^^^'^/^ 
are still left-right hybrids but each has higher amplitude on 
its representative surface. Thus, we have conveniently sep- 
arated the degenerate surface states in a way that they natu- 
rally evolve into the single surface state case in the limit of 
infinite film thickness. It should be stressed here that maxi- 
mizing I (Sk) surf I is only the means of choosing a reasonable 
approximation to the single-surface state. However, we an- 
alyze the resulting spin-polarization field by calculating the 
spin expectation value of ^|^^^'^ over all layers, (Sk) = 

l^^max,L|^|^max,L^^ 

In Fig. [5] we present our results on the spin-polarization 
field, represented by arrows, for an 1 1 -layer W(OOl) film. The 
figure is focused on the area of reciprocal space around the M- 
point, i.e., partly outside the first Brillouin zone, to show the 
contours which correspond to the surface states SI and S2. 
The starting point of each arrow corresponds to a k-point on 
the Fermi surface, the direction of the arrow to the direction 
of (Sk), and the color code to its magnitude. The out-plane 
component of spin polarization is found to vanish zero. 

The spin-polarization field of surface state S 1 (two concen- 
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trie rings around M) reminds one of the "pure" Rashba states 
at the Au(lll) surface, in the sense that (Sk) is in-plane and 
almost perpendicular to k — ko, where ko is the ring cen- 
ter (here the M-point; in the Rashba model for Au(lll) the 
r point). More precisely we observe that (Sk) is perpen- 
dicular to the group velocity, i.e., (Sk) is along the Fermi- 
surface tangent, as expected from the Rashba model. As for 
the spin-polarization field of the state S2, first we should re- 
mind the reader that S2 merges with the bulk- state continuum 
(see Fig. [2]), excluded from the plot, which leads to the fact 
that the shown contour in Fig. [5] ends seemingly abruptly. 
Second, S2 also has a partner with opposite direction of the 
spin-polarization field, which lies higher in energy (see the 
band- structure in Fig.|4]). (Sk) in S2 is also in-plane and per- 
pendicular to the group velocity. 

Finally, it should be noted that the magnitude of the 
spin-polarization fields, | ^ (^k (r) | | ^k (r) ) | cannot reach the 
maximal value ^ , except perhaps in the pure Rashba model, 
since an interaction with other states at other energies at the 
same k via the spin-orbit operator will always be present 
and will reduce the value. To view this from a different 
perspective, the states ^k(^) ^^^^^ linear combination 
^max(^j.j = ak^k (^) + /^k^k (r) result in a non-collinear 
spin-density ^^k(r)^^''^k(r) that can be brought in a diago- 
nal form only in an r-dependent reference frame; but in order 
to achieve | (Sk) | = ^, this reference frame would have to be 
independent of r. 




FIG. 6: (Color online) Value of b^{s) (top) and T~^{s) (in units 
of ps~^/at%) for s on the unit sphere for the case of a 10-layer 
W(OOl) film with a W adatom as scattering defect. The highest spin- 
relaxation values are found for s out-of-plane (here taken as the z- 
axis). 



IV. ANISOTROPY OF THE ELLIOTT- YAFET 
PARAMETER AND OF THE SPIN-FLIP SCATTERING 
RATE OFF ADATOMS 

We now examine the anisotropy of the spin-relaxation rate 
with respect to the polarization direction of the injected spin 
population (i.e. the SQA s) relative to the crystallographic 
axes. For bulk systems we have already discussed the effect in 
Ref. 14, pointing out that the reduced symmetry of thin films 
will clearly play a role in this phenomenon. To summarize 
the origin of the effect, even though the spin-orbit operator 
L • S is independent of the SQA, its matrix elements are not. 
Particularly relevant here are the matrix elements of the spin- 
flip SOC {LSy^, \{L+S- + ^-5+), that are responsible for 
the spin-mixing parameter 6^ and for the spin-flip transitions 
in general. Due to this dependence on the SQA, and consid- 
ering the crystallographic symmetry of the W(OOl) film, we 
expect at least three inequivalent directions of s for which the 
spin-relaxation rate T^^ will become extremal: perpendicu- 
lar to the film, in-plane in the [100] direction and in-plane in 
the [110] direction. Comparing the values of T^^{s) in all 
directions we obtain the definition of the anisotropy 

A[T-^] = [md.xT-\s) - mmT-\s)]/mmT-\s) (12) 

s s s 

which is a somewhat different quantity compared to A[b'^] 
introduced in Ref. 14, where we had b'^{s) in the place of 
T^^{s) implying the Elliott approximation. We wish to point 
out here that not only the spin-flip rate but also the spin- 



conserving rate depends on the SQA, and that part of the 
anisotropy comes from the SQA dependence of 6^ in the 
Bloch states, while another part comes from the spin-flip scat- 
tering off the impurity potential. The total scattering rate (spin 
flip plus spin conserving) defined by Eqs. ^ and ([TO]) is inde- 
pendent of s. 

In Fig. [6] we show in a color code the value of b'^{s) and 
T^^{s) as a function of the direction s on the unit sphere for 
a 10-layer W(OOl) film with a W adatom as scatterer. The 
maximum spin-relaxation rate of T^^ = 11.63 ps~^/at% 
is obtained for s out of plane; the minimum, T^^ = 7.91 
ps~^/at%, for s along the [1 10] axis. This yields an anisotropy 
of ^[T]~^] = 47%. The map of b'^{s) also shows a clear max- 
imum for s out of plane with 6^ = 0.294 while in the [110] 
direction it has a value of 6^ = 0.215 giving an anisotropy of 
A[b'^] =37%. Evidently there is no complete quantitative cor- 
relation between b'^{s) and T^^{s), since the Elliott approxi- 
mation is too crude in this case, but qualitatively the correla- 
tion is obvious. It should be noted that the in-plane variance 
of either b'^{s) or T^^{s) is small (on the order of 2%), ow- 
ing to the high symmetry of the fourfold crystallographic axis. 
The anisotropy here for the thin films has a significant higher 
value compare to the anisotropy of 6^ of about 6% in bulk W, 
as we have found previously, due to the reduced symmetry 
and to a great extent due to the surface states 

We can analyze the effect further by examining the Fermi- 
surface-resolved 6^ and T^~^(k) [defined in Eq. ([t])]. Fig.]?] 
shows these quantities for s || [001] (left) and 5 || [110] (right). 
Both quantities [6^ and T^~^(k)] show a symmetry compati- 
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s II [001] s II [110] 




FIG. 7: (Color online) Distribution of the spin-mixing parame- 
ter bl- (top) and spin-relaxation rate T^'^{]<.) (bottom) shown in a 
color-code (in units of ps~^/at%) on the Fermi surface of a 10- 
layer W(OOl) film. For the calculation of T^"^ (k), scattering off W 
adatoms was considered. For the calculations shown in the left pan- 
els the spin-quantization axis was taken in the [001] direction (i.e., 
perpendicular to the film), for the calculations shown in the right 
panels it was taken in the [110] direction (in-plane). The distribu- 
tions on the Fermi surface are clearly compatible with the specific 
crystal symmetry operations that either leave s unchanged or result 
in s —s. Especially on the surface states the values of 6^ and 
Tj"^ (k) are highest in the regions where the spin-polarization fields 
(Fig.|5]) are perpendicular to s. 



V. REMARKS ON THE RELAXATION MECHANISM IN 
THE RASHBA STATES 

Finally we wish to discuss the physical mechanism of spin 
relaxation in the presence of Rashba states. In semiconduc- 
tors or semiconductor heterostructures it is naturally assumed 
that the spin-polarization fields contribute to spin relaxation 
via the Dyakonov-Perel mechanism.^ ^ There, an electron oc- 
cupies a state at k that is split in energy by only a little ac- 
cording to an effective Hamiltonian H = —^fl\^-a, thus 
in a semiclassical picture the electron spin precesses around 
the spin-orbit field ft\^ with a Larmor frequency |r2k|, losing 
memory of its original direction before being scattered away 
(within this effective model the spin polarization field is par- 
allel or antiparallel to the spin-orbit field ft\^). However, for 
this to happen it must be assumed that t he el ectron wavepacket 
has an energy spread larger than | fik | This is possible in 
semiconductors where \ is usually small^ (of the order of 
ImeV or less, depending on temperature, doping concentra- 
tion, etc.), as usually oc k with k very close to the con- 
duction band minimum at F. As opposed to this, in metal 
surfaces the Rashba splitting A Ef of surface states at can 
be large, of the order of 100 meV, e.g. A^;^ ^ 200meV in W 
(Fig.0), 30meV in Cu(lll),^ and ISOmeV in Au(lll).I^It 
is unlikely that a coherent wavepacket excited by e.g. an in- 
jected spin current or by microwave radiation in CESR should 
have such a large energy spread, activating the Dyakonov- 
Perel mechanism, except perhaps if prepared very precisely 
by an experiment targeting exactly this. Therefore, we con- 
sider the Dyakonov-Perel mechanism not applicable in the 
case of the Rashba surface states of most metals, but we can- 
not exclude it for special cases e.g. lighter metals with signifi- 
cantly weaker spin-orbit coupling or favorable band structure 
(e.g. in Ref. 42 we see that Ag(lll) shows a splitting of only 
= 2meV due to the shallow surface state). 



VI. SUMMARY AND CONCLUSIONS 



ble to the intersection of the symmetry operations leaving the 
Fermi surface invariant and the operations leaving the SQA 
invariant plus the inversion. It is interesting to see that even 
though the system parameters are outside the prerequisites of 
the Elliott approximation (6^ is rather large and the scatter- 
ing off a transition metal adatom cannot be considered weak), 
still a correlation between the "highs and lows" of 6^ and 
T]~^(k) is visible in the color code. Also, that the surface 
state regions show high values for both quantities when the 
spin-polarization fields (see Sec. Ill and Fig. |5]) are perpen- 
dicular to the SQA and lower values when they are parallel 
or antiparallel to the SQA. In this sense, and since there will 
be a part of the Rashba states with spin-polarization field per- 
pendicular to the SQA no matter what the choice of the SQA 
is (compare with Fig. [5]), the surface states are clearly "spin 
drains" for the system due to the spin polarization fields that 
are related to the Rashba effect; the worst case is s || [001], 
when all spin-polarization fields are perpendicular to s. 



In summary, we have investigated spin-relaxation physics 
of W(OOl) ultrathin films from first principles. We observe 
that the Elliott- Yafet parameter exhibits an even-odd oscilla- 
tion with respect to the number of layers of the films, which 
stems from an even-odd oscillation in the surface electronic 
structure of the films at and near the Fermi energy. The oscil- 
lation is then inherited by the spin-relaxation rate that depends 
on the Elliott- Yafet parameter. 

We have further identified the Rashba character and spin- 
polarization fields of the surface states and discussed how they 
contribute to the anisotropy of the spin-relaxation rate with re- 
spect to the relative orientation between the spin-quantization 
axis and the crystallographic directions. The anisotropy val- 
ues are much higher compared to those in bulk W. 

We believe that our findings are not only particular to 
W(OOl) free standing films but are more general at least for 
transition metals in the bcc structure and even when sand- 
wiched between insulators. We base this speculation on three 
considerations: First, we performed calculations (not pre- 
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sented in the present paper) of the ElHott-Yafet parameter for 
Mo(OOl) films in the bcc structure and found basically the 
same oscillatory effect. Second, the existence of surface states 
is, in general, closely connected to the crystal structure; all 
bcc transition metals will show a dip in the density of states 
at the center of the (i-band and in all cases the breaking of 
translational symmetry at the (001) surface will produce sur- 
face states of the character found here within the gap of the 
surface-projected bulk band-structure. Third, the mechanism 
just stated is expected to produce interface states in the case 
that the film is in contact with an insulator, even if the details 
of the interface band structure can be more complicated in this 
case. 

As an outlook we believe that it is worthwhile to investi- 



gate these effects in a broader family of ultrathin films. Such 
investigations are in progress and will be reported in a future 
publication. 
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